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Abstract 

Basic chiral QCD parameters (the pion decay constant, the quark and gluon 
condensates, the dynamically generated quark mass, etc) as well as the vac¬ 
uum energy density (up to the sign, by definition, the bag constant) have 
been calculated from first principles within a recently proposed zero modes 
enhancement (ZME) model of the true QCD vacuum. Our unique input data 
was chosen to be the pion decay constant in the chiral limit as given by the 
chiral perturbation theory at the hadronic level (CHPTh). In order to ana¬ 
lyze our numerical results we set a scale by two different ways. In both cases 
we obtain almost the same numerical results for all chiral QCD parameters. 
Phenomenological estimates of these quantites as well as vacuum energy den¬ 
sity are in good agreement with our numerical results. Complementing them 
by the numerical value of the instanton contribution to the vacuum energy 
density, we predict new, more realistic values for the vacuum energy density, 

the bag constant and the gluon condensate. 

PACS numbers: 11.30 Rd, 12.38.-t, 12.38 Lg and 13.20 Cz. 
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I. INTRODUCTION 


Let us begin the second part of our paper with the discussion of one of the most interesting 
feature of (dynamical chiral symmetry breaking) DCSB. As was underlined in the first part 
(hereafter referred to as I), there are only five independent quantities by means of which all 
other chiral QCD parameters should be calculated. For the sake of convenience, let us write 
down them together. They are: 


F 2 - —h 2 -- 1 
r CA ~ 8n 2 K o z o 


[ Zo . zB 2 (z 0 ,z) 

'o {zg 2 (z) + B 2 (zo, z)} ’ 


m d = k o {z o B 2 (z o ,0)} 1/2 , 


( 1 . 2 ) 


(qq)o = -^J 0 dzzB(z 0 ,z), (1.3) 

zq 

e q = ~^ k o z o 2 f dz z {In z [zg 2 (z) + B 2 (z 0 ,z)] -2zg(z) + 2}, (1.4) 

o 

e g = -\k A 0 z^ 2 x [l8 ln(l + y) - ^z 2 ln(l + —) - • (1-5) 

TT L O 2 Zq 2 J 

Recall that g(z) and B 2 (zq, z) are explicitly given by (1.14) and (1.15) of I. It is instructive 
also along with them to write down definition (3.13) of I for DCSB scale, namely 

A csBq = 2m d . (1.6) 

So these final expressions which should be used to calculate chiral QCD parameters within 
our approach depend only on two independent quantities, namely: mass scale parameter ko 
and the constant of integration of dynamical quark SD equation of motion z 0 . However, 
it follows from (1.2) that information on the parameter Zo should be extracted again from 
m d and the initial mass scale parameter k o itself, which characterizes the region where 
confinement, DCSB and other nonperturbative effects begin to play a dominant role (see 
below). So the second indepent parameter zq is reduced to the pair of the mass scale 
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parameters ko and Despite the fact that in onr consideration the initial mass scale 
parameter /j (characterizing the scale of nonperturbative effects) has been introduced by 
” hand”, such a transformation of pair of independent parameters ko and zo into the pair of ko 
and rrifi is also a direct manifestation of the phenomenon of the ” dimensional transmutation” 
[1], This phenomenon occurs whenever a massless theory acquires masses dynamically and 
it is a general feature of spontaneous symmetry breaking in held theories. 

Let us emphasize once more that it generally characterizes our approach, in order to 
calculate numerically all chiral QCD quantities (considered here and others), that one needs 
only two independent (free) parameters, both having significant and clear physical sense. 
The above mentioned dynamically generated quark mass m^, playing a role of the integra¬ 
tion constant of the corresponding equation of motion (the quark SD equation) instead of 
because of the above mentioned ” dimensional transmutation” and mass scale parameter ko, 
responsible for a scale at which important nonperturbative effects begin to play a dominant 
role. Our calculation scheme is self-consistent because we calculate n = 5 independent phys¬ 
ical quantities by means of m = 2 free parameters having clear physical sense, so condition 
of self-consistensy n > m is satisfied. The general behaviour of all of our parameters as 
given by the relations (1.1-1.5) are shown in Figs. 1-5. 

Our approach makes it possible to calculate all chiral QCD parameters (the ones consid¬ 
ered here plus others) at any requested combination of and k 0 , but in order to analyse 
numerical results it is necessury to set a scale at which it should be done. We set a scale 
by two, at first sight, different ways but leading (see below) to almost the same numerical 
results within our calculation scheme. 

Evidently, to set a scale in each case makes it possible to determine only one of the 
two free parameters in our calculations. In order to determine the second one we use the 
chiral value of the pion decay constant obtained by the chiral perturbation theory at the 
hadronic level (CHPTh) in Ref. 2, namely F° = (88.3 ± 1.1) MeV. This value is chosen 
as an input data in our numerical investigation of chiral QCD. The pion decay constant 
is a good experimental number since it is a directly measurable quantity as opposed, for 
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example, to the quark condensate. For this reason we may reliably estimate the deviation 
of the chiral values of physical quantities, which can not be directly measured, from their 
” experimental”, phenomenologically determined values, when the chiral value of the pion 
decay constant is approximated by the experimental value. 

In the above mentioned CHPTh ( or equvivalently the effective held theory) there is a 
low energy constant B, determined by (qq) 0 = — F 2 B and measures the vacuum expectation 
value of the scalar densities in the chiral limit. It is just this constant that determines the 
meson mass expansion in the general case. Indeed, in leading order (in powers of quark 
masses and e 2 ) from CHPTh, one has [3, 4] 


Ml + = 

('™°u + rn° d )B 

(1.7) 

m 2 k+ = 

(■ m° u + m° s )B 

(1.8) 

M 2 k o = 

(m° d + m° s )B. 

(1.9) 


Calculating (independently) the constant B, one then will be able to correctly estimate 
current quark masses rn () u , m d and rn° s by using the experimental values of meson masses [5] 
in (1.7-1.9). 

II. ANALYSIS OF THE NUMERICAL DATA AT A SCALE OF DCSB AT THE 

QUARK LEVEL 

Let us begin by recalling that there exists a natural scale within our approach to DCSB. 
Indeed, at the fundamental quark level the chiral symmetry is spontaneously broken at a 
scale A-csBq defined by (1.6). Therefore it makes sense to analyse our numerical data at a 
scale where DCSB at the fundamental quark level occurs. To this end, it is necessary only 
to simply identify mass scale parameter ko with this scale A csBq, he. to put 

ko = A csBq = 2 m,d- (2.1) 

In other words, we will analyse our numerical results at a scale responsible for DCSB at the 
fundamental quark level. Evidently, this uniquely determines the constant of integration of 
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the quark SD equation. Indeed, from (2.1) and on account of (1.2), then it immediately 
follows that this constant is equal to Zq = 1.34805. From the pion decay constant in the 
chiral limit, chosen as input data, and on account of this value for zq, from (1.1) it yields the 
numerical value for /e 0 . This means that all physical parameters considered in our paper are 
uniquely determined. As it was mentioned above, it will be instructive to explicitly display 
our numerical results when the chiral value of the pion decay constant is approximated by 
the experimental value advocated in Refs. 6 and 7, namely F° = 92.42 MeV , as well as by 
the standard value F° = 93.3 MeV. Results of our calculations are displayed in Table 1. 

Let us make a few concluding remarks. To set a scale by the way described in this section 
has the advantage that it is based on the exact definition (1.6) for a scale of DCSB at which 
analysis of the numerical data must be done. In general, it is not obvious that this scale 
A csBq and scale A c , at which quark confinement occurs, should be of the same order of 
magnitude. Moreover, the information about A c is hidden within this scheme of calculation. 
In order to reveal the raison d’etre for A c and its relation to A csBq, let us set a scale in the 
way described in the next section. 

III. ANALYSIS OF THE NUMERICAL DATA AT THE CONFINEMENT SCALE 

In our approach there exists only one scale, denoted as /j or k 0 (separating, in general, 
the nonperturbative phase from the perturbative one), which is responsible for all the non¬ 
perturbative effects in QCD at large distances. If there is a close relation between quark 
confinement and DCSB (and we believe that this is so) then the scale of DCSB at the funda¬ 
mental quark level (1.6) and the confinement scale A c should be, at least, of the same order 
of magnitude. In other words, in our approach A c should be very close to A csBq- This is in 
agreement with Monte Carlo simulations on the lattice which show that the deconfinement 
phase transition and the chiral symmetry restoring phase transition occur approximately at 
the same critical temparature [8], hereby confirming the close intrinsic link between these 
nonperturbative phenomena. 
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Unfortunately, neither the exact value of 'rn c i nor of ko is known. For this reason, let us 
first reasonably assume that the dynamically generated quark masses, in any case, 

300 < rn d < 400 (MeV), (3.1) 

should hold but otherwise they remain arbitrary. We believe that this interval covers all 
possible realistic values used for and obtained in various numerical calculations. The second 
independent parameter ko should be varied in the region of 1 GeV - the characteristic scale 
of low energy QCD. Varying independently these pairs of parameters rn d and ko numerically, 
one can calculate all chiral QCD parameters with the above derived formulae. 

From the value of the pion decay constant in the chiral limit, as well as from the range 
selected first for rn d (3.1) and on account of (1.1) and (1.2), it follows that the momentum k 0 
always should satisfy the upper and lower boundary value conditions, namely 691.32 < k 0 < 
742.68 {MeV). The vacuum energy density contributions of the nonperturbative gluons 
(1.5) changes its sign in the range selected for rn d (3.1) and in this interval for ko- Therefore 
it becomes positive and this should not be so because of the normalization condition (we 
normalize perturbative vacuum to zero). It is easy to show that this is result of that that 
the lower bound chosen for the dynamical generated quark mass in (4.1) is too low. Indeed, 
the vacuum energy density (1.5) vanishes at the critical point z c q = 1.45076 (see Fig. 6). 
Then from (1.2) calculated at this point, it follows that 

ko < 2.26 rnd- (3.2) 

Using this ineqaulity in additional, the vacuum energy density (1.5) will always be negative 
as it should be and it will become zero only at critical values determined as k^ = 2 . 26 md- 
From the chosen interval for md (3.1) and the obtained interval for ko, however, it follows 
that the ratio between the corresponding lower bounds ko/rud = 691.32/300 = 2.3044 does 
not satisfy the above obtained inequality (3.2), while this ratio for the corresponding upper 
bounds ko/rrid = 742.68/400 = 1.8567 satisfies it. This explicitly shows that the lower bound 
for m d in (3.1) was incorrectly chosen. The exact lower bound for rrid can be found from the 
ko r as 742.68 = 2.26 md, and (3.1) becomes 
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328.62 < m d < 400 (MeV). 


(3.3) 


In the range determined by (3.3) and in the above obtained interval for ko, the vacuum 
energy density (1.5) will be always negative because any combination (ratio) of ko and rn d 
from these intervals will satisfy inequality (3.2). But this is not the whole story yet. A new 
lower bound for m d leads to a new lower bound for ko as well. Indeed, combine now this 
new lower bound (3.3) with the chiral value of the pion decay constant one obtains a new 
lower bound for ko as well. 

As noted above, ko is regarded as a momentum which separates the nonperturbative 
phase (region) from the perturbative one. In the region obtained for ko the nonperturbative 
effects, such as quark confinement and DCSB, begin to play a dominant role. It is a region 
determining a scale at which confinement occurs. From now on let us call this scale for ko 
a confinement scale (in the chiral limit) and denote it A c . So the final numerical values for 
the confinement scale are as follows 

707 < A c < 742.68 (Me V). (3.4) 

In intervals determined by (3.3) and (3.4) the vacuum energy density e g will be always 
negative (see Fig. 7). 

It is worth noting that any value for A c from interval (3.4) is possible but not any 
combination of A c from interval (3.4) and rn d from interval (3.3) will automatically satisfy 
the value of the pion decay constant. Therefore it is necessary to adjust values of m d from 
(3.3) for chosen value of A c from interval (3.4) and vice versa (see Fig. 8). This means that 
m d is in close relationship with A c . Moreover, completing the above mentioned procedure, 
one finds that A c is nearly the double of the generated quark mass m d , i. e. A c « 2 m d . This 
confirms that A c and A csBq defined by (1.6) are nearly the same indeed. In the previous 
calculation scheme the adjusting procedure was automatically fulfilled because of the exact 
relation (2.1). Thus there is an intimate relationship between A csBq and A c on the one hand 
and the double generated quark mass m d on the other hand. 
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The interval (3.4) for possible values of A c along with the new range for rri c i (3.3) will 
uniquely determine numerically the upper and lower bounds for all other chiral QCD param¬ 
eters considered here. Like in the previous case, our numerical results are shown in Table 
2 (calculation scheme B), where the shorthand (0|G 2 |0) stands for the gluon condensate 
(0|^G“ I/ Cr“j0). Our numerical bounds for the vacuum energy density e need additional 
remarks. We note that the bounds for e is not the sum of bounds for e q and e g . The upper 
and lower bounds for e q are achieved at the upper and lower bounds for rn ( j (A c ) while for 
e g they are achieved at the lower and upper bounds of rn ( j (A c ). 

Let us now prove the relation A c « 2 rrid- We have already learnt that correct values of ko 
belong to the interval for A c (3.4). Then identifying with A c in (3.2), one finally obtains 

A = ±(—1 + -——) < 0.13, (3.5) 

A CSBq 

where the positive sign corresponds to A c > A csBq and the negative one is valid when 
A c < AcsBq■ In the derivation of this relation we used definition (1.6). 

Finally it is worth underlining once more that besides good numerical results obtained in 
this section, we have established the existence of realistic lower bound for the dynamically 
generated quark masses. In each calculated case their numerical values are shown in Table 
2. Thus one concludes that the vacuum energy density due to the nonperturbative gluons is 
sensitive to the lower bound for m^. Another important result is that we have clearly shown 
that the confinement scale A c and DCSB scale A csBq are nearly the same indeed. 

IV. CONCLUSIONS AND DISCUSSION 

Let us briefly compare our numerical results obtained from first principles with phe¬ 
nomenologically estimated values of the physical parameters considered here. An estimate 
of the quark condensate in Refs. 9 and 10, 

(qq.)1 /3 = ~ 


(225 ± 25) MeV 


(4.1) 



is in good agreement with our values. It is worth noting here that QCD sum rules give 
usually the numerical values of physical quantities, in particular the quark condensate, 
approximately within an accuracy of (10-20)% (see, for example Ref. 11). 

Our values for the current quark masses are also in good agreement with recent estimates 
from hadron mass splittings [12] 

ml = (5.1 ±0.9) MeV, 
ml = (9.0 ± 1.6) MeV, 

m° s = (161 ± 28) MeV (4.2) 

and QCD sum rules [13] 

m° u = (5.6 ± 1.1) MeV, 
m° d = (9.9 ± 1.1) MeV, 

m° s = (199 ± 33) MeV, (4.3) 

see also reviews [14]. 

Here it is worth noting that from our numerical results (Tables 1 and 2) it follows that the 
constituent quark mass m q should differ little from m^. Apparently, the difference between 
them is of the order of a few per cent only of the displayed values of m^. So without 
making a big mistake even for light quarks, it is possible to simply use rrid instead of rn q . 
Doing so one comes to the conclusion that the CHPTh value of the pion decay constant 
and the constituent quark model (CQM) with the value for the constituent quark mass 
rn q = 362 MeV advocated by Quigg [15] are nearly in one-to-one correspondence within our 
calculation scheme (see Table 1). Moreover, from our numerical results (Tables 1 and 2) one 
can conclude that the dominant contributions to the values of all chiral QCD parameters as 
well as the vacuum energy density come from large distances, while the contributions from 
the short and intermediate distances can only be treated as small perturbative corrections. 

The phenomenological analysis of the QCD sum rules [10] for the numerical value of the 
gluon condensate implies 
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(4.4) 


<0|^G^Gyo>~ 0.012 GeV 4 , 

and using then (2.5) of I, one obtains the vacuum energy density as 

e ~ -0.003375 GeV 4 . (4.5) 

In the random instanton liquid model (RILM) [16] of the QCD vacuum, for a dilute ensemble, 
one has 

e = ~ x 1.0 f m~ 4 ~ -0.003411 GeV 4 . (4.6) 

The estimate of the gluon condensate within the QCD sum rules approach can be changed 
within a factor of two [10]. We trust our numerical results for the vacuum energy density 
much more than those of the gluon condensate. The former was obtained on the basis of 
the completely nonperturbative ZME model of the vacuum of QCD while the latter was 
obtained on account of the perturbative solution for the CS-GML /3-function [10]. Let us 
also emphasize the one important fact that our calculation of the vacuum energy density 
is a calculation from first principles while in the RILM [16] the parameters characterizing 
vacuum, the instanton size p o = 1/3 fm and the ” average separation” R = 1.0 fm were 
chosen to precisely reproduce traditional (phenomenologically estimated from the QCD sum 
rules) values of quark and gluon condensates, respectively. 

We reproduce values (4.4-4.6), which are due to the instanton-type fluctuations only, 
especially well when the pion decay constant in the chiral limit was approximated by its 
experimental value. Moreover, our numerical results clearly show that the contribution to 
the vacuum energy density of the confining quarks with dynamically generated masses e q is 
approximately equal to e g , that is of the nonperturbative gluons. It is well known that in 
the chiral limit (massless quarks) tunneling is totally suppressed, i.e. the contribution of the 
instanton-type fluctuations to the vacuum energy density vanishes. It will be restored again 
in the presence of DCSB [17-19]. Thus, in principle, in the chiral limit and in the presence 
of DCSB, the total vacuum energy density should be the sum (as minimum) of these three 
quantities, i.e. 
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G — £l + £g + Nf€q, 


(4.7) 


where 6/ describes the contribution of the instanton component to the vacuum energy density. 
We introduce also the explicit dependence on the number of different quark flavors Nf since 
e q itself is the contribution of a single confining quark. Of course, this should be valid for 
the non-chiral case as well. The distinction will be in the concrete values of each component, 
apart from, maybe, e g . 

Let us now run a risk and make a few quantitative predictions. Indeed, not going into the 
details of the instanton physics (well described in Ref. 18) and in agreement with the authors 
of Ref. 10, it is worth assuming, for simplicity’s sake, that the light and heavy quarks match 
smoothly. This allows one to choose for the instanton component of the vacuum energy 
density ej the average value between (4.5) and (4.6), i.e. namely ej ~ —0.0034 GeR 4 . Then 
our predictions for more realistic values of the total vacuum energy density e t (for Nf light 
confining quarks with dynamically generated masses) and the corresponding values of the 
gluon condensate are listed in Tables 3 and 4, 5 for both calculation schemes A and B, 
respectively. 

It is worth reproducing explicitly some interesting particular values of the total vacuum 
energy density and the corresponding values of the gluon condensate. Thus for a pure 
gluodynamics (Nf = 0) one has 

e t ~ —0.005 GeR 4 , 

(0|^G^G;j0)~0.0177GeV, (4.8) 

and 

- 0.00661 <e t < -0.003837 (GeR 4 ), 

0.0136 < (0|f G“„G“J0>< 0.0235 (GeR 4 ). (4.9) 

Here and below the numbers correspond to the approximation of the pion decay constant in 
the chiral limit by its standard value. 
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For the more realistic case Nf = 2 one obtains 

e t ~ -0.008 GeV 4 , 

( o \^ g ; m g ; jq ) ~ 0.0233 Get/ 4 , (4.10) 

and 

- 0.00933 <e t < -0.00724 (GeV 4 ), 

0.0256 < (0|^G^G“J0)< 0.0331 (GeV 4 ). (4.11) 

There exist already phenomenological estimates of the gluon condensate [20] as well as 
lattice calculations of the vacuum energy density [21] pointing out that the above mentioned 
standard values (4.4) and (4.5-4.6) are too small. Our numerical predictions are in agreement 
with these estimates though we think that their numbers for gluon condensate [20] are too 
big. At the same time, it becomes quite clear why the standard values are so relatively 
small, because they are due to the instanton component of the vacuum only. 

In the above mentioned RILM [16], light quarks can propagate over large distances in 
the QCD vacuum by simply jumping from one instanton to the next. Within our model (see 
I and Ref. 20) propagation of all quarks is determined by the corresponding SD equations 
(clue to the ZME effect) so that they always remain off mass-shell. Thus we need no picture 
of jumping quarks. Contrast to the RILM, we think that the main role of the instanton-like 
fluctuations is precisely to prevent quarks and gluons from the freely propagation in the 
vacuum of QCD. Running against instanton-like fluctuations quarks undergo difficulties in 
their propagation in the QCD vacuum which, in principle, is a very complicated inhomoge- 
nious medium. At some critical value of the instantons density the free propagation of the 
virtual quarks, apparently, become impossible, so they already never annihilate again with 
each other. Obviously, this is equivalent to the creation of the quark-antiquark pairs from 
the vacuum. From this moment nontrivial rearrangement of the vacuum can start. The 
above mentioned critical value can be reached when ej ~ e g + e q , i.e. when at least one 
sort of quark flavors is presented in the QCD true vacuum. On one hand, this is supported 
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by our numerical results for e = e g + e q . On the other hand, the numerical values for ej, 
as given by (4.5) or (4.6), also confirms this. In the realistic (nonchiral) case the instanton 
part, along with other contributions, may substantilly differ from those shown in Tables 1, 
2, 3, 4 and 5. 

The bag constant is defined as the difference between the energy density of the pertur¬ 
bative and the nonperturbative QCD vacuums. We normalize the perturbative vacuum to 
zero (I). So in our notations the bag constant becomes 

B = -e t (4.12) 

(Not to be mixed with the CHPTh constant (1.7-1.9)). Our predictions for this quantity 
are also shown in Tables 3 and 4, 5 for each calculation scheme A and B, respectively. In 
fact, our values for the bag constant overestimate the initial MIT bag [23] volume energy by 
one order of magnitude. Nevertheless, we think that the introduction of this constant into 
physics was a main achievement of the bag model. 

As in previous case, let us explicitly reproduce some interesting concrete values of the 
bag constant. For a pure gluodynamics (Nf = 0) it is: 

B ~ 0.005 GeV 4 ~ (266 MeV)) 4 ~ 0.651 GeV/fm 3 . (4-13) 

For the more realistic case Nf = 2 the bag constant becomes 

B ~ 0.008 GeV 4 ~ (300 MeV) 4 ~ 1 GeV/fm 3 . (4.14) 

For simplicity’s sake we reproduced its value obtained within the calculation scheme A only. 
It has been noticed in [24] that noybody knows yet how big the bag constant might be, but 
generally it is thought it is about 1 GeV/fm 3 . The predicted value for Nf = 2 is in fair 
agreement with expectation. 
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TABLES 


TABLE I. Calculation scheme A 


F° 

TV 

88.3 

92.42 

93.3 

MeV 

A CSBq 

724.274 

758.067 

765.284 

MeV 

rn d 

362.137 

379.0335 

382.642 

MeV 

(QQ) o 

(-208.56) 3 

(-218.29) 3 

(-220.36) 3 

MeV 3 

e i 

-0.0012 

-0.00143 

-0.0015 

GeV A 

e 9 

-0.0013 

-0.00157 

-0.0016 

GeV A 

e 

-0.0025 

-0.0030 

-0.0031 

GeV 4 

(0|^G^G^|0) 

0.009 

0.0106 

0.011 

GeV 4 

m° u 

6.65 

6.36 

6.30 

MeV 

m° d 

10.08 

9.63 

9.54 

MeV 


202.85 

193.75 

191.94 

MeV 
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TABLE II. Calculation scheme B 


F° = 88.3 

F° = 92.42 

F° = 93.3 

707 < A c < 742.68 

737.9 < A c < 768.4 

744.4 < A c < 773.86 

328.62 < m d < 400 

340 < m d < 400 

342.416 < rn d < 400 

(-210.34) 3 < (qq) Q < (-206.9) 3 

(-219.3) 3 < (qq) Q < (-216.34) 3 

(-221.2) 3 < (qq) 0 < (-218.33) 3 

-0.00135 <e q < -0.00096 

-0.0016 < e q < -0.00128 

-0.0017 < e q < -0.00136 

-0.0024 < e g < -0.00045 

-0.00226 < e g < -0.00044 

-0.00221 <e g < -0.000437 

-0.00336 < e < -0.0018 

-0.00354 < e < -0.002 

-0.00356 < e < -0.0021 

0.0064 < (0|G 2 |0) < 0.0128 

0.007 < (0|G 2 |0) < 0.0192 

0.00746 < (0|G 2 |0) < 0.0199 

6.48 < < 6.81 

6.27 < ml < 6.53 

6.22 < ml < 6.47 

9.83 < m° d < 10.33 

9.5 < m° d < 9.89 

9.43 < m° d < 9.81 

197.67 < m° s < 207.7 

191 < m° s < 199 

189.76 < m° s < 197.34 
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TABLE III. Calculation scheme A. Predictions 


F° 

7T 

92.42 

93.3 

MeV 

e* = e/ + £g + A T /e q 

-0.00497 - A/0.00143 

-0.005 - A/0.0015 

GeV 4 

(0\^G^G- u \0) 

0.01767 +A/0.00508 

0.01777 + A/0.00533 

GeV 4 

B 

0.00497 + A/0.00143 

0.005 +A/0.0015 

GeV 4 
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TABLE IV. Calculation scheme B. Predictions 


F° = 92.42 

-0.00566 - iV/0.00128 < e t < -0.00384 - N f 0.0016 
0.0136 + iV/0.00568 < (0|G 2 |0) < 0.020 + iV/0.0045 
0.00384 + JV/0.0016 < B < 0.00566 + iV/0.00128 


20 







TABLE V. Calculation scheme B. Predictions 


F° = 93.3 

-0.00661 - 7V/0.00136 < e t < -0.003837 - TV/0.0017 
0.0136 + 1V/0.006 < (0|G 2 |0) < 0.0235 + 1V/0.0048 
0.003837 + 1V/0.0017 < B < 0.00661 + 1V/0.00136 
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FIGURES 

FIG. 1. The pion decay constant Fqa as a function of ko, drawn only for the most reasonable 
region, selected first for the dynamically generated quark masses (3.1). The obtained interval for 
ko is also explicitly shown (see Section 3 below). 

FIG. 2. The quark condensate as a function of ko, drawn only for the most resonable region, 
selected first for the dynamically generated quark masses (3.1). 

FIG. 3. The vacuun energy density due to confining quarks with dynamically generated masses, 
as a function of ko, drawn only for the most resonable region, selected first for the dynamically 
generated quark masses (3.1). 

FIG. 4. The vacuum energy density due to the nonperturbative gluons as a function of ko, drawn 
only for the most resonable region, selected first for the dynamically generated quark masses (3.1). 
The obtained interval for ko is also explicitly shown (see Section 3 below). 

FIG. 5. The vacuun energy density e as a function of ko, drawn only for the most resonable 
region, selected first for the dynamically generated quark masses (3.1). 

FIG. 6. The vacuum energy density due to the nonperturbative gluons contributions (1.5) as 
a function of zo- 

FIG. 7. The vacuum energy density due to the nonperturbative gluon contributions e g as a 
function of ko- A c is the confinement scale (3.4). A new interval for m c j (3.3) is also shown. A 
similar figure can be drawn for the case when the pion decay constant is approximated by the 
experimental value. 

FIG. 8. The pion decay constant Fqa as a function of ko- A c is the confinement scale (3.4). A 
new interval for rn c ] (3.3) is also shown. A similar figure can be drawn for the case when the pion 
decay constant is approximated by the experimental value. 
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